EPJ manuscript No. 

(will be inserted by the editor) 



Directed motion of domain walls in biaxial ferromagnets under 
the influence of periodic external magnetic fields 
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Abstract. Directed motion of domain walls (DWs) in a classical biaxial ferromagnet placed under the influ- 
ence of periodic unbiased external magnetic fields is investigated. Using the symmetry approach developed 
in this article the necessary conditions for the directed DW motion are found. This motion turns out to 
be possible if the magnetic field is applied along the easiest axis. The symmetry approach prohibits the 
directed DW motion if the magnetic field is applied along any of the hard axes. With the help of the soliton 
perturbation theory and numerical simulations, the average DW velocity as a function of different system 
parameters such as damping constant, amplitude, and frequency of the external field, is computed. 

PACS. 05.45.Yv Solitons - 75.10. Hk Classical spin models - 75.78. Fg Dynamics of domain structures 



1 Introduction 

One-dimensional ferromagnetic models are currently of 
considerable experimental and theoretical interest [TJ[2]. 
A large portion of this interest has been directed towards 
the domain wall (DW) response to an ac magnetic field. 
One of the important problems here is the development 
of different ways for obtaining a net DW drift under the 
influence of unbiased perturbations. 

The ratchet effect [3j[H[5j|6] has been shown to be an ef- 
ficient tool to control the motion of particles and particle- 
like excitations. The mechanism of this effect is based on 
the breaking of all symmetries that connect two solutions 
with specular velocities [7]. For topological solitons this 
phenomenon has been investigated both theoretically [HI 
' l!flll0lllllll2[[T3l and experimentally [T41IT51H6] in continu- 
ous and discrete Klein-Gordon-type systems. It has been 
shown 10,11 that a biharmonic external field, consisting 
of a sinusoidal signal and its even overtone can yield a 
directed soliton motion. Similar biharmonic external field 
has been used in |17) to control the dynamics of an indi- 
vidual spin. 

It should be mentioned that spatial asymmetry can be 
used for controlling the domain wall motion. This can be 
achieved by creating a sawtooth-like asymmetric pattern 
on the magnetic film [18,19 . Experimental observation of 
this phenomenon was reported in |20j . Asymmetric pin- 
ning potential that consists of triangular holes has been 
proposed and experimentally implemented in [21]. Obser- 
vations have shown that this asymmetry favours certain 
direction of the domain wall propagation. 



On the other hand, since the work of Schlomann [33J 
[23] the problem of a DW drift under the influence of an os- 
cillating magnetic field, polarised either in the plane con- 
taining the easy axis [2"4"ll2"51l2"rj] or in the plane perpen- 
dicular to it 22. 23. 123], has been studied in the literature. 
Despite a certain number of papers devoted to this prob- 
lem, an interesting and important question arises: what 
are the necessary conditions which one has to impose on 
the unbiased external periodic magnetic field, such that a 
unidirectional DW motion will arise as a result? Also, we 
would like to point out the need of a unifying approach 
that would join together different ways of driving a uni- 
directional DW motion. In this paper, we show that the 
symmetry approach [TUlHTj should be a perfect tool for 
this task. 

Thus, the aim of this work is to investigate in detail 
the possibility of the unidirectional motion of magnetic 
topological solitons (domain walls). In particular, we for- 
mulate the necessary conditions which have to be imposed 
on the external unbiased magnetic field in such a way that 
this motion will take place. 

The paper is organised as follows. In the next section, 
we describe the equations of motion for the biaxial fer- 
romagnet. In Section 3, the symmetries of the Landau- 
Lifshitz (LL) equation are discussed. The average domain 
wall velocity is computed analytically in Section 4. The 
numerical solution of the LL equation is given in Section 
5. Conclusions and a final discussion are presented in the 
last section. 
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2 Model and equations of motion 

The dynamics of the one-dimensional chain of classical 
spins in the continuum limit is described by the well- 
known Landau-Lifshitz (LL) equation 



9 t S = -[Sx (9*S + JS)] + ef(S,t), 



(1) 



where S(x,t) = (S x , S y , S z ) is a three-component di- 
mensionless magnetisation vector. Without loss of gen- 
erality, we assume the following normalisation condition: 



o ,, 



-S* = 1. The matrix J = diag( J x , J y , J z ) contains 



the information about the anisotropy constants (/3i = 
Jx — J y , k = Jz — Jy), so that the total energy of the 
magnet is given by 



E= l - 

2 



-Hoc 



[(d x S) 2 - SJS] dx = 



+oo 



£[S{x,t)]dx. (2) 



Note that if /3i < 0, 03 > \@i\, the OZ axis is the easiest 
axis, so that we have an easy axis ferromagnet with XY 
being the anisotropic hard plane. Here £(S) is the energy 
density function. The perturbative term ef contains the 
external magnetic field and the phenomenological Landau- 
Gilbert damping 



ef=-[SxH(t)]+A[SxS t ]. 



(3) 



Here the periodic external magnetic field H(t) = H(i + 
T) has zero mean value (H(t)) t = and A is a damping 
constant. For the most of magnetic materials A ~ 0.008 -f- 
0.01 [27|. 

In this paper, we are interested in computing the aver- 
age DW velocity as a function of system parameters. But 
before embarking on this task, it is necessary to investi- 
gate the symmetry properties of the LL equation. 



3 Symmetries of the Landau-Lifshitz equation 

According to the previous work [TTlfTU] , we state that the 
necessary condition for the occurrence of the directed DW 
motion is the breaking of all the symmetries that relate 
two solitons with the same 

topological charge and with specular velocities: 



S S(x, t; v) = S(x, t; —v). 



(4) 



The unperturbed (e = 0) DW solution of the LL equa- 
tion is well known [28]. It is a topological soliton S' ' (x, t; <j>) 
= (cos <f> sechz, sine/) sechz, Qtanhz), where z = £(x + 
xo) ~ Q(t=£(x - x - vt), £ = y/(3 3 ~ Pi cos 2 cj>, C = 
/3i sin 4> cos 4>- Here Q is the topological charge of the DW 
soliton with Q = 1 corresponding to the kink (soliton) 
and Q = — 1 to the antikink (antisoliton) solution. The 
azimuthal angle <fi describes the direction of the projection 
of the magnetisation vector S on the XY plane, and thus 
the handedness or polarity of the DW changes depending 
on the interval to which the value of <j> belongs: < < 7r 



or — 7r < (f> < 0. The DW velocity is defined by the value 
of (f>, moreover, v(<p) — v(4> + tt) and t>(— <ft) = — v ((/)). 

Taking into account the properties of the unperturbed 
soliton and assuming that the perturbation is weak enough 
not to distort the soliton shape, one can define the soliton 
velocity and the center of mass as 



dX c 

dt 



1 

E 



+ 00 



x^£[S(x,t)]dx. (5) 



It is easy to see that there exist only two types of symme- 
tries that can relate two arbitrary solutions with opposite 
velocities and the same topological charge. These oper- 
ations must include either a space reflection and a time 
shift or, vice versa, a time reflection and a space shift: 

o i a : x } x -r x , t ^ t H- T a , o a ^ »_> a , 

S z -t-S z ; a = x,y, (6) 

S2a ■ X — > X + x' , t — >• — t + t' , S a — >• — S a ', 

a = x,y, (7) 

where x' and t' are arbitrary constants and t q = nT/2, 
n = 0,1,2. For the sake of clarity, let us briefly discuss 
the symmetries S\ x and Si y . The symmetry Si x acts on 
the unperturbed solution by turning a DW S^°\x,t;<j)) 
into a solution S^°'(x,t;TT — <j>), while the symmetry S\ y 
turns a solution with (f> into a solution with —<f>. In these 
cases, v(ir — (f>) = —v{4>) and v(— <f>) = —v(<j>), respec- 
tively. Therefore, the abovementioned symmetries connect 
two DW solutions with opposite velocities, while the other 
DW properties, in particular, the topological charge Q or 
the width £ _1 , remain unchanged. Note that although the 
symmetry Si y changes the polarity of the DW, we still 
consider these two solutions as the same because the po- 
larity is a local characteristic of a DW, in contrast to the 
topological charge. 

Application of the perturbation ([3]) can or cannot de- 
stroy the above symmetries. Note that when the magnetic 
field is applied, the energy density @ must be comple- 
mented with the term — (S,H(i)). The symmetries iSi a 
are present if there exists such t that the following equal- 
ities are satisfied for the magnetic field H(t): 

H a (t + Ta ) = ~H a (t), Hp{t + T a ) = Hf>(t), 

H z (t + Ta )=-H z (t); a = x,y, f3 ^ a . (8) 

One should stress that for each symmetry Si a there is a 
corresponding value of r a . 

The symmetries S^ a are always violated in the pres- 
ence of dissipation (A ^ 0), however, if A = 0, they are 
present if there exists such t' that the following equalities 
take place: 

H a (-t + t') = -H a {t), Hfi(-t + if) = H {t), 

H z (-t + t')=H z {t); a = x,y,(3^a. (9) 

Thus, in a general case A ^= 0, in order to obtain the 
directed soliton motion one has to apply a magnetic field 
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for which in both (for a — x and a — y) the sets of the 
equalities (JSJ at least one equation in each set does not 
hold. In the dissipationless case, one has to violate both 
the sets © and ©. 

Consider now the oscillating magnetic field directed 
along one of the coordinate axes. If H(i)||OZ, in order 
to break the symmetries, S\ a one has to choose the re- 
spective component H z (t) in such a way that it satis- 
fies the inequality H z (t + T/2) ^ —H z {t). In the cases 
H(t)||OX and H(t)||OY, one of the symmetries S\ a will 
always be present. Indeed, if H(t) = (H x (t),0, 0) T , the 
symmetry S± x is present if there exists such t x that H x (t+ 
t x ) = —H x (t). Obviously, this equation holds if t x — 
T/2. The symmetry S\ v is present if there exists such 
t v that the equality H x (t + r y ) — H x (t) holds. Since 
the external magnetic field is periodic, this condition is 
automatically fulfilled for r y = T . Similarly, if we con- 
sider H(t) = (0,iJ y (t),0) T , the symmetry S\ y is present 
if H y (t + T/2) = —H y (t), but for the presence of Si x it 
is sufficient to guarantee the periodicity of the function 
H y (t). Thus, it is not possible to obtain a directed soli- 
ton motion by applying magnetic field only along OX or 
OY axis, at least, for arbitrary small perturbation. On 
the other hand, it is possible to obtain a directed motion 
in the case H(t)||OZ, since one can violate the equality 
H z (t + T/2) = —H z (t) by various choices of the magnetic 
field, for instance, by choosing it in the following bihar- 
monic form : 

H(i) =e z H(t) =e z [Hicos(u)t) + H 2 cos(mut + 9)}, (10) 

with m = 2, 3, . . . . If H% ^ and m is even, the above- 
mentioned equality is always violated, while for odd m's 
it is always satisfied. In the dissipationless case A = 0, an- 
other set of symmetries, namely S201, must be broken. For 
a given choice of the magnetic field direction, this situation 
occurs if there does not exist such t' that H z {—t + t') = 
H z (t). For the function (fTU|) this means that the symme- 
tries S2 a are violated if m is even and H 2 ^ 0, 9 7^ 0, ±7r. 

4 Computation of the average DW velocity 
using the perturbation theory 

The perturbation theory in the first order implies that 
the perturbation ([3]) is too small to distort the soliton 
shape and it influences only the temporal evolution of the 
following soliton parameters: its center of mass X(t) and 
the azimuthal angle 4>(t). The equations that describe the 
parameter evolution are obtained in accordance with the 
papers [29,30 : 



where f a (a — x, y, z) are the components of the per- 
turbation vector ([3]), and the soliton parameters £(t) — 

VPs ~ A cos2 </H^) an d C{t) = Pi sin (/>(£) cos 4>(t) now de- 
pend on time. 

Consider the case of the magnetic field (jlOl) directed 
along the OZ axis. In this case, the equations for the time 
evolution of the soliton parameters X(t) and 4>(t) have the 
following form: 



-f=H(t) + rsm2<j), T = 

at 

dX _ /3isin2(/. 

dt ~ 2^3-^icos 2 ^' 



A/3i 
2 : 



(13) 
(14) 



If one assumes that the soliton shape changes only in- 
significantly due to the perturbation, the vibrations of the 
azimuthal angle 4>(t) around its equilibrium position ap- 
pear to be small. Next, we assume the amplitudes H\2 
to be small parameters. Then the oscillating solution of 
the equation flT3"|) should be sought in the form 4>(t) = 
±tt/2 + fli0i(t) + Hl4> 2 {t) + Hf</) 3 (t) + Q(Hf). The ini- 
tial value ±7r/2 comes from the fact that the Bloch wall 
configuration is energetically most favourable. After rep- 
resenting the external magnetic field in the form H(t) = 
Hih(t) = Hi[coscut + H2/H1 cos (2uit + 8)} and substitut- 
ing the expansion for cj)(t) into equation (|13l) . one finds 
an approximate expression for 4> which contains only the 
(t>i(t) and 4>3(t) terms. Then it remains only to substitute 
the expression for 4>(t) into equation (|14[) and to average it 
over one oscillation period. In the expansion of the r.h.s of 
the equation (fT4")l into the Taylor series with respect to the 
parameter H\ , we have limited ourselves with the term of 
the order O(Hf). As a result, the following expression for 
the average DW velocity is obtained: 

(v)~-QA(H li2 ,u,r)sm(d-6o), 

3/3? H 2 Hl 



A[Hi fl ,u,r) 



16/3| /2 (4r 2 



2 )VT2 



8q = 2a,rcta,n(2T/uj) — arctan(_T/w). 



(15) 



We would like to stress that this expression is the same 
for the initial angles <f> = 7r/2 and <\> — —ir/2. 

In the case of odd m's, we obtain (v) = 0. The ex- 
pression (I15p clearly confirms the validity of the symme- 
try approach. The average soliton velocity becomes zero if 
H 2 = and this signals the restoration of the symmetries 
Si a . In the dissipationless limit (A — >• 0), 9 — >• and thus 
(v) ex sin 9. In this limit, the symmetries S201 are restored 
if H(t) = H(—t). The average DW velocity becomes zero 
at the values 9 = 0, ±tt, so that they are precisely those 
values at which the function H{t) is symmetric. 



d(f> 

dt 


e r +co 


coshz 


dz, 


dX 
~dT ' 


" £ 2f 


/•+00 

J — 00 


£ 2 coshz 


f-(z) = 


= f x sin 4> - 


fy COS 


0, 



Qfz 



( ' 5 Numerical simulations 

dz, In order to verify the symmetry approach developed in 

Section [21 the initial LL equation has been discretized in 

(12) the spatial dimension with the step h = 0.05 and the 
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resulting system of coupled ordinary differential equations 
has been integrated numerically using the fourth order 
Runge-Kutta method. It is convenient to solve numerically 
the LL equation in the following form: 



-ftS = 



1 



1 + A 2 
H, = d*S + JS 



[S x H*] 
H(t), 



A 



[S x [S x H,]] , 



(16) 



which is equivalent to equation ([1} with the perturbation 
([3]). The validity of this method has been checked by mon- 
itoring the energy conservation in the purely Hamiltonian 
case A = and H(t) = 0. 

It should be emphasized that in order to compute 
the mean DW velocity, one has to average over the set 
of initial conditions: phase — 7r < <p < 7f, initial time 
< to < T, and time t. In the numerical simulations, 
we consider only the dissipative case A > 0, therefore we 
are interested in attractor(s) that correspond to moving 
DWs. If the perturbation Q is small, then the phase space 
of the LL equation have to consist of the basin(s) of at- 
traction of periodic attractor(s) [limit cycle(s)] that are 
locked to the frequency w of the external magnetic field. 
Breaking the respective symmetries should manifest itself 
in desymmetrization of the basins of attraction that cor- 
respond to DWs moving with opposite velocities. Below 
we demonstrate that actually in the case of broken sym- 
metries there exists only one attractor that corresponds 
to the directed DW motion. In this case, it is sufficient 
to compute only the average velocity on the attractor: 
(v) -^t^oo [X(t + T) - X(t)]/T, where X(t) is the DW 
center of mass. 

First, we consider the case H(i)||OZ and the expres- 
sion for the Z-component of the magnetic field given by 
equation (fTO"]) . The time evolution of the DW center com- 
puted as X max = max a;e (_ 00!+00 )f [S (*,*)] (shown in Fig- 
ure [TJ clearly demonstrates the validity of the symmetry 
approach. Curve 1 in this figure corresponds to the case of 
the single harmonic drive (H2 — 0), where no directed DW 
motion is seen. Next, no directed DW motion is observed 
in the case of mixing two odd harmonics (m = 3) as shown 
by curve 2. Curves 3-5 illustrate the evolution of the DW 
center for the case of two mixed harmonics with m = 2. It 
is easy to notice that on the time scale t > A -1 , the sys- 
tem settles on a periodic attractor (the periodicity can be 
observed from the insets in Figure [TJ) that corresponds to 
the motion in the direction, defined by the phase shift 9. 
The simulations have been performed for different initial 
values of and initial times to, and in all the cases the sys- 
tem settles on the same attractor (compare, for example, 
curves 4 and 5). 

The dependence (v) (9) of the average DW velocity on 
the phase shift between the harmonics, 9, is shown in Fig- 
ure [5J It appears to have a sinusoidal shape, as predicted 
by the perturbation theory result (TTSl) . Another demon- 
stration of the validity of the symmetry approach is the 
behavior of the points where (v)(6) — 0. In Figure [5] along 
with the data for A = 0.01, the data for A = 0.1 and 
A = 0.3 have been plotted as well. Although these values of 
damping do not correspond to realistic values for magnets, 




Fig. 1. Time evolution of the kink center X max under action 
of the external magnetic field (|10|l for J x — 0.2, J y — 0.5, 
J z — 1, A = 0.01, id — 0.2. Curve 1 corresponds to the case of 
single harmonic drive (Hi = 0.1, H2 = 0), curve 2 corresponds 
to H 1 = H2 — 0.1, m = 3, 9 — —2. Curves 3-5 correspond to 
Hi=H 2 = 0.1, m = 2, 6 = -1 (curve 3), 9 = 2.5 (curves 4,5). 
The last two curves show evolution for different initial times: 
to = T/A (curve 4) and to = (curve 5). Insets (a)-(c) give 
details of curves 1,4, and 5, respectively. 




Fig. 2. Dependence (v) = (v)(9) for A = 0.01 (©), A = 0.1 (o, 
scaled by factor 1.3), and A = 0.3 (♦, scaled by factor 1.85) 
computed numerically from equation (|16[) . The perturbation 
theory result (|15|l for A = 0.01 is shown by solid line. Other 
parameters are chosen the same as for curves 3-5 in Figure [Tj 
The inset shows details in the vicinity of 9 = 0. 



these results are very instructive for the illustration of the 
restoration of the symmetries S2 a - Indeed, when A — > 0, 
the values of 9, at which the DW velocity becomes zero, 
gradually shift to 0,±7r. As shown already in Section [3J 
both the symmetries S2a are restored if H z (—t) = H z (t). 
This happens if 9 — 0,±7r [see equation (|T0l) ]. We would 
like also to point out good correspondence between the 
perturbation theory results (shown by the solid line) given 
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by equation (fT5|) and the results of direct integration of 
the LL equation. 

Next, we consider the case of the magnetic field di- 
rected along the OX axis [H(t) = e x H(t)], but which has 
the same functional dependence (ITU1) . According to the 
symmetry approach, this periodic drive does not yield the 
directed DW motion. The numerical simulations of this 
case are illustrated by Figure[3] The basins of attraction of 
two DW solutions that have opposite polarities and oppo- 
site velocities appear to be symmetric with respect to the 
value <j) = 0. Indeed, if the initial value of the azimuthal 
angle is positive (see curves 1,3, and 5), the dynamics of 
the system settles on the attractor with the positive DW 
velocity (v = 0.00031), while for negative values of 4> (see 
curves 2,4 and 6) it tends to the solution with the DW 
velocity of opposite sign (v — —0.00031). In the case of 




Fig. 3. Time evolution of the kink center X ma x under ac- 
tion of the external magnetic field (|10p for H\ = H2 = 0.05, 
6 = 1, m = 2 with initial azimuthal angle (j> — 7r/10 (curve 
1), <f> = — 7r/10 (curve 2), <f> — 7r/4 (curve 3), 4> — — 7r/4 (curve 
4), (j> — 37r/4 (curve 5), and 4> = — 37r/4 (curve 6). The rest of 
parameters is as in the case of Figure [T] The inset shows distri- 
bution of the magnetisation vector components (S x , S y ,S z ) T 
of DWs that correspond to curves 3 (black) and 4 (red) of the 
main figure. 



H(t)||OY, the same scenario is observed. 



6 Discussion and conclusions 

In this paper, the symmetry approach for the analysis of 
the unidirectional domain wall motion has been developed. 
The main objective was to demonstrate that a proper 
choice of the oscillating unbiased magnetic field can yield 
a net directed motion of the domain wall. With the help of 
the symmetry approach, we have obtained the necessary 
conditions to be imposed on the magnetic field in order 
to obtain the unidirectional motion. When the magnetic 



field is applied along a certain coordinate axis, this mo- 
tion turns out to be possible only if this axis coincides with 
the easy axis, say OZ. The symmetry arguments prohibit 
the directed DW motion if the magnetic field is applied 
along any of the hard axes. The necessary condition for 
the directed DW motion in the presence of dissipation is 
given by H z {t) 7^ —H z (t+T/2). Next, we demonstrate the 
cubic dependence (v) ex H\H2 of the average DW veloc- 
ity on the magnetic field amplitude. These results appear 
to be similar to those for topological solitons in the ac 
driven sine-Gordon (SG) equation [ID]. Very recently, [3T] 
a rigorous proof of the universality of these results for a 
wide range of nonlinear systems driven by the biharmonic 
signals of the type (fTU)) has been obtained. 

However there are certain differences in the directed 
soliton motion in the LL and SG cases. Due to the fact that 
the LL equation is three-component, while SG is scalar, 
there are more ways to apply the external field to the 
system, namely along any of three coordinate axes. But 
only in the case of magnetic field applied along the easy 
axis, the directed motion is possible. Here it should be 
stressed that we are interested in the average net motion 
which is independent from the the local properties of the 
DW such as handedness (polarity) . The application of any 
ac signal along non-easy axes drives DWs with opposite 
polarities into opposite directions. The basins of respective 
attractors of the LL equation are symmetric with respect 
to the 7r-shift of the initial azimuthal angle. Therefore, a 
weak noise which is inevitable in realistic systems will lead 
to exploration of the whole phase space and eventually to 
zero net motion. 

Another difference with respect to the SG case caused 
by the multicomponentness of the LL equation is a wider 
range of possibilities to drive a DW by external oscil- 
lating fields. Let us briefly discuss the case when two of 
the magnetic field components are nonzero. Consider first 
the case of magnetic field polarised in the easy plane XY 
[H(t) = {H x (t),H y (t),0) T }: H x (t) = H^ } cos(m x ojt + x ), 

H y (t) = Hq cos(m y ujt + 6y). Here m x ,m y — 1,2,3,..., 
and they do not have a common divisor. Note that such a 
way to control the DW motion has been suggested in [21] 
for the particular case of m x = m y = 1. The symmetry 

S\ x requires the simultaneous fulfilment of the equalities 

H x (t+T/2) = -H x (t) &ndH y (t+T/2) = H y {t). Similarly, 

the S\y symmetry is present if both the equalities H x (t 4- 
772) = H x (t) and H y (t + T/2) = -H v (t) hold. If both m x 
and m y are odd, we have H a (t + T/2) = —H a (t) 7^ H a (t). 
Therefore none of the sets of equations © can be satisfied 
and thus both symmetries Si a are broken. If m x is even 

H x (t) ^ -H x {t), 
therefore 



■Hy(t), 



and m y is odd, we have H x (t + 772) 

thus S lx is broken. But, H y (t + T/2) = 

Si y is present. Similarly, if m y is even and m x is odd, the 

symmetry Si x is satisfied and Si y is broken. Therefore if 

one of m Xl y is even and another one is odd, we expect no 

directed motion, whereas this motion must occur if both 

Tn x , v are odd. 

Another way to drive unidirectionally a domain wall 
is to apply the oscillating magnetic field polarised in the 
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plane that contains the easy axis |26j : 

H y (t) = H^ v} cos(m y ujt + 6 y ), H z (t) = H^ z) cos(m z ujt + 



0; 



1, 2,3,.. .. Alternatively, one can consider 



26 



the field polarised in the XZ plane. If both m v and m z 
are odd, it is impossible to break simultaneously the set 
of equalities (JHJ) and (0) . In this case, the breaking of S\ x 
takes place because H z (t + T/2) = —H z (t) but H y (t + 

T/2) = —H y (t) ^ H y (t). However, the symmetry S\ y 
is still present. If m y is even and m z is odd, we obtain 25 
H z (t+T/2l=-H z (t) but H y (t+T/2) = H y (t) ^ -H y (t), 
therefore S\ x is present while Si y is broken. If m y is odd 

and m z is even both the symmetries Si x . y are broken be- 
cause H z (t + T/2) = H z (t) ^ -H z (t). Therefore this is 
the only way to obtain the directed DW motion with the 
help of magnetic field polarised in the YZ plane. 

Finally, we would like to outline the future directions 
of applications of the symmetry approach. It is of interest 
to consider topological magnetic excitations in two- and 
three-dimensional systems, where alongside the directed 
translational motion a unidirectional rotation can take 
place as well, as shown previously for particles [32]. An- 
other question is how to apply the symmetry approach to 
the problem of DW directed motion in more complicated 
magnetic systems such as antiferromagnets, ferrites, mag- 
netoelastic systems, and others. In this direction, some 
progress has already being accomplished in the papers [531 
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